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■"^.^ , Context. Bisectors of strong, single spectral lines, usually the Fei6252 line, have traditionally been used to examine the velocity fields in 
r^^ stellar atmospheres. This requires high S/N often achieved by summing many individual spectra, 
tin Aims. We investigate whether bisectors derived from cross-correlation functions (CCF) of single-exposure spectra can be used to provide 
<»>>, information on stellar atmospheres, and whether they can be used to discriminate between radial velocity changes caused by planets, magnetic 
5_j activity and oscillations. 

^ Methods. Using a sample of bright stars observed with the HARPS spectrograph, we examine the shapes of the bisectors of individual strong 

C^ spectral lines in summed spectra, comparing with similar studies in the literature. Moreover, we examine four different quantitative CCF 

kJ bisector measures for correlations with radial velocity and stellar parameters. 

• ,-H Results. We show that CCF bisector measures can be used for quantitative analysis, employing both the absolute values and the variations. 

ji^ > From absolute values, log g and absolute magnitude can be approximated, and from the correlations with radial velocity one can distinguish 

?H ' between magnetic activity, oscillations and orbiting planets as the probable cause of radial velocity variations. We confirm that different 

. . . I isolated spectral lines show different bisector shapes, even between lines of the same element, calling for caution in trying to derive global 

stellar properties from the bisector of a CCF. For the active star HR 1362 we suggest from the bisector shape an extra photospheric heating 

caused by the chromosphere of several hundred degrees. We confirm the fill-in of spectral lines of the Sun taken on the daylight sky caused by 

Rayleigh-Brillouin and aerosol scattering, and we show for the first time that the fill-in has an asymmetric component. 

Key words. Instrumentation: spectrographs - Techniques: radial velocities - Line: profiles - Stars: atmospheres - planetary systems - Stars: 
activity 

1. Introduction the granulation and its cause, the convection, is harder to in- 
terpret in terms of line profile variations. These variations are 

The analysis of high-resolution stellar spectra can reveal a ^f^^^ described in terms of the line bisector. The physical in- 

wealth of information about the physical conditions in the terpretation of the line asymmetries giving rise to the different 

atmosphere. Immediate information on effective temperature, ^-^^^^^^ ^j^^p^^ j^^^ ^^^^ discussed by "o^ ("200?) from an ob- 

gravity and elemental abundances comes from the analysis of servatio nal point of vie w, while theoretical modelling has been 

Ime depths and equivalent widths (EW), while information on ^^^^ bv lAst)lund et alJ (Eoool) with good results for the Sun. 

T , „ ~~ ~ Z TT T-. 11 1 4.J m<a The bisector has diagnostic power for a number of stellar 

Send offprint requests to: T. H. Dall, e-mail: tdall@eso . org ° ^ 

* Based on observations collected at the La Silla Observatory, ESO Parameters: The effective temperature and luminosity can be 

(Chile), with the HARPS spectrograph at the ESO 3.6m telescope read from the shape of the bisector (the classical "C" shape) 

(ESO Programme 075.D-0760), and on observations obtained from and from the height of the blue-most point on the bisector, 

the ESO/ST-ECF Science Archive Facility (ESO Programmes 073.D- respectively (.Grav.^OOSLl . Changes in the bisector shape may 

0578, 073.C-0784 and 074.D-0380). have several different causes: Variations at a very low level 
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can be introduced by reflected light from an orbiting planet 
llHatzes et al.lll998l) . More importantly, the bisector is heavily 
afifected by photospheric magnetic fields, and as such it has 
been used to distinguish radial velocity (RV) changes caused 
by planets from changes in the p hotospheric velocity fields in- 
duced by magnetic activity (e.g. lOueloz et alJl200lh . The bi- 
sector has been used to infer the presence of unseen stellar 
companions, both physically connected compan ions as well 
as unrelated objects along the line of sight (e.g. ISantos et alJ 
[2002; Torres et al. 2005, respectively). Classical stellar oscil- 
lations also introduce spectral line asymmetries that will be 
visible in the bisector, but rather than the traditional bisector, 
other mea sures have often b een used fo r oscillation studies (e.g. 
iBaldrv & Bedding 2000; D aU & Frandserj2002l) . 

High precision single-line bisector studies require very high 
S/N ratios, which for bright stars can be achieved with a single 
exposure. For fainter stars, or if the instrument does not offer 
sufficiently high S/N in one exposure, the high S/N is often 
achieved by summation of several tens to hundreds of individ- 
ual spectra. If the summation is done with spectra taken over 
a long period of time, it may introduce errors due to devel- 
opments in the atmosphere (e.g. oscillations or activity) and 
due to errors on the wavelength calibration of individual spec- 
tra. Even small wavelength offsets will have an impact on the 
summed spectral line and its bisector. This is also true for the 
case where bisectors are calculated for many different lines in a 
single exposure and then averaged, since the absolute positions 
of the individual bisectors are generally not known with very 
high precision. 

The recent advent of a new generation of ultra-stable high- 
resolution spectrographs has caused a revolution in the detec- 
tion of extra-solar planets in the past few years. Since the first 
detection of an extra-solar planet around a solar-like star (51 
Peg; Mayor & Oueloz 1995), the art of detecting planets by 
measuring RV displacements has been steadily refined. The 
methods used are the self-calibration with an iodine cell and 
the simultaneous ThAr calibration method. While the former 
is self-calibrated due to the passing of the stellar light through 
iodine vapor, the second relies on the recording of a ThAr cal- 
ibration spectrum alongside the stellar spectrum. The simulta- 
neous ThAr method has the advantage that the stellar spectrum 
is undisturbed, while the iodine-cell introduces hundreds of I2 
absorption lines in the spectrum. On the other hand, the ThAr 
method imposes much stricter requirements on the mechani- 
cal stability of the spectrograph. The cuiTent state-of-the-art 
instrument for accurate RV measurements using the simulta- 
neous ThAr m ethod is beyo nd doubt the ultra-stable spectro- 
graph HARPS (! Mavoretal.L.200 3). instaUed at the 3.6m tele- 
scope at the La Silla site of the European Southern Observatory, 
which can ach ieve precision better than 1 ms ' per exposure 
JRupprecht et al...2004) . The high precision of HARPS comes 
from the high internal mechanical, temperature and pressure 
stability, combined with the use of the simultaneous record- 
ing of a ThAr calibration lamp spectrum alongside the stel- 
lar spectrum. In the reduction process, the calibrated spectrum 
is cross-correlated with an appropriate stellar mask, matched 
to the spectral type of the target star. The resulting cross- 
correlation function (CCF) can be thought of as a "mean" 



spectral line, and as such, it has inherently high S/N. Iodine 
cell spectra have also been used for CCF bisector analysis by 
iMartinez Fiorenzano et al.l 112005^ after proper removal of the 
iodine lines. However, this implies an extra reduction step, and 
in what follows we will only consider CCF's obtained with the 
simultaneous ThAr method. Due to the high spectral stability, 
calibration errors are negligible, and errors due to changes with 
time are eliminated since the CCF is constructed from one sin- 
gle exposure, allowing better time resolution to study variabil- 
ity phenomena. Furthermore, in contrast to the classical bisec- 
tors, an absolute velocity position of the HARPS CCF bisector 
is known to very high precision. On the other hand, combining 
many different spectral lines, spanning a wide range in wave- 
length and physical conditions may wash out any physical in- 
formation, and it does not discriminate between blended and 
unblended lines. Moreover, the particular mask used in con- 
structing the CCF will have an effect on the shape of the bisec- 
tor. 

In this paper we will investigate the bisectors measured 
with HARPS on late-type stars. Especially, we will investigate 
whether the CCF can be used as a meaningful probe of stellar 
atmospheric physics instead of the classical bisectors of indi- 
vidual spectral lines normally used. In order to do this we com- 
pare the classical bisector found from the sum of tens to hun- 
dreds of individual spectra with the CCF bisector for a small 
grid of solar-like stars, both active and non-active. 



2. Observations and data reduction 

The HARPS pipeline (Data Reduction Software; DRS) gives 
the reduced ID spectrum plus a file containing the computed 
CCF for all 72 orders, in addition to the mean CCF, with the 
computed RV included as a header keyword. 

In order to compare the CCF bisector with the classical bi- 
sector, we have combined a large number of ID spectra for 
each target star (Table Q in order to reach as high S/N as 
possible. The wavelength calibration and stability of HARPS 
are good enough that negligible errors are introduced due to 
the adding of spectra taken on different dates. The observa- 
tions were obtained during long asteroseismology runs, dur- 
ing planet searches and during technical tests of the guiding 
system. The solar spectra were taken on the daytime sky for 
calibration purposes during one of the asteroseismology pro- 
grammes (Kjeldsen et al., in preparation). 

Each individual spectrum was first shifted to the rest wave- 
length using the calculated RV and the IRAF task dopcor, and 
then combined using scombine. 



3. Bisector definitions and tools 

The ge neral definition of the bisector has been given bv lGravi 
( 1988). In short, one finds the midpoint of the line for a number 
of intensity positions inside the line. In practice, one may calcu- 
late the midpoint between all the points in the blue (red) wing 
and corresponding interpolated points in the red (blue) wing. 
For our analysis we interpolate both the red and blue wing in 
turn, which results in a better sampling of the bisector. 
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Table 1. The target stars. A^sp is the number of ID spectra 
combined. The S/N was estimated from several continuum 
regions in the 6200-6600 A range of the combined spectra. 
Atmospheric parameters are from Santos et al. (2005) (HR 98, 
6585, 7665), Allende Priet o et alJ (|2004 (Sun, HR 509, 2261, 
5933, 7602), Dal l et ah i,2005ii) (HR 1362) and this work (HR 
1326. see Sect. lO. 



HR Name 



Spec. 
Type 



Kr 



S/N 



[K] 



logg [Fe/H] 





Sun 


G2V 


628 


1200 


5777 


4.44 


0.00 


98 


/?Hyi 


G2IV 


2766 


1700 


5837 


4.00 


-0.08 


509 


rCet 


G8V 


55 


1500 


5328 


4.62 


-0.52 


1326 


a Hor 


Kim 


137 


800 


4670 


2.75 


+0.02 


1362 


EKEri 


G8III 


9 


1300 


5240 


3.55 


+0.09 


2261 


a Men 


G6V 


36 


1000 


5473 


4.51 


+0.03 


5933 


ySer 


F6IV 


160 


1300 


6246 


4.30 


-0.15 


6585 


fi Ara 


G3V 


275 


1100 


5806 


4.28 


+0.32 


7602 


/3Aql 


G8IV 


135 


1200 


5106 


3.54 


-0.19 


7665 


5Pav 


G7IV 


76 


700 


5614 


4.29 


+0.36 



Once the bisector has been found, several parameters can be 
extracted regarding the shape. Tradit ionally, the so-call ed 7%- 
span or t he velocity s pan is used cToner & Gravlll988l) . More 
recently, lOravl (120051) found that the height of the blue-most 
point on the bisector is a powerful luminosity indicator for the 
cooler spectral types, while the temperature dependence in- 
volves the full shape of the bisector and must be probed by 
co mparing with t abulated values. 

lOueloz et alJ ( l2001.) pioneered the use of bisectors of CCFs. 
To distinguish between planet- and activity-induced RV varia- 
tions, they calculated a velocity span, which they note is es- 
sentially the inverse slope of the bisector and analogous to the 
velocity span definition used for classical single-line bisector 
work. We will reserve the term velocity span for classical bi- 
sectors, and we will adopt the definition of Queloz et al. of the 
bisector inverse slope: 



BIS = Vt - V},, 



(1) 



where Vt is the mean bisector velocity in the region between 
10% and 40% of the line depth, and Vb is the mean bisector ve- 
locity between 55% and 90% of the line depth. This definition 
is also similar, but n ot identical, to the velocity span defined by 
|PovichetalJ(l200ll) . 

While the BIS is a way to quantify the bisector, it is a purely 
arbitrary measurement. In this paper we will compare the BIS 
with three other bisector measurements, that are defined as fol- 
lows: (1) The bisector slope b/,, which is defined as the inverse 
slope from a linear fit to the part of the bisector between 25% 
and 80% of the line depth. The reason for this choice, is that the 
CCF bisectors often appear very straight in the central parts, as 
opposed to the characteristic "C" shape of a classical single- 
line bisector (2) The curvature, defined as 



Ch = (V'3 - V2) - (V2 - Vi), 



(2) 



which follows IPovich et alJ (1200 ih . except that we define the 
mean velocities in segments including 20-30% (vi), 40-55% 




200 m/s 



Fig. 1. The Fei /1 6252 bisectors of all program stars. Crosses 
show the calculated points of the bisector with eiTor bars on the 
velocity. 

{v2), and 75-100% (V3) of the line depth. We do not use the ve- 
locity displacement Vh of Povich et al. because it is essentially 
just an average bisector velocity position, relative to the posi- 
tion of the line core A^. However, this normalization is purely 
arbitrary. A better choice would be to measure it relative to the 
line center i.e. the RV, which in general does not correspond to 
the line core position. Instead we employ (3) the bisector bot- 
tom Vbot, which is simply the average of the bottom four points 
on the bisector minus the RV. This is indeed as arbitrary a con- 
struction as any, and can only be justified by considering that 
we know the absolute position of the line (the RV) to very high 
precision. 

Here we must define what we mean by the absolute posi- 
tion of the line: While the true (in an absolute sense) value of 
the RV is not known, we can measure absolute positions rel- 
ative to the instrument zeropoint, which is given with 1 ms"^ 
precision in this case, defined by the long-term stability of the 
HARPS spectrograph. This precision depends on the S/N, so 
for individual lines it will be somewhat lower than for the CCF. 

4. Bisectors of single lines 

The bisectors of the Fei A 6252 line are plotted in Fig. Q] Error 
bars are estima ted from the measured S/N in the spectra, using 
the formula of Gray ( 1988). This is the only line used by Gravl 
(2005), and comparing the stars we have in common (J3 Aql, 
T Cet and the Sun) we find very good overall agreement. Hence 
HARPS can — not surprisingly — be used as a classical high- 
resolution spectrograph for line-profile analysis, although it is 
necessary to sum several individual HARPS spectra to reach 
the high S/N. 

Although in an overall sense, the bisector shapes can be un- 
derstood in terms of simple temperature and luminosity effects. 
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Table 2. All individual spectral lines investigated. 



Line 




;r[eV] 


Notes 


Fei 


6151.62 


2.17 


Fig.lA.ll 


Fei 


6252.56 


2.40 


Fig.E] (1) 


Cat 


6499.65 


2.52 


(2) 


Nil 


6643.63 


1.67 


Fig.O 


Fei 


6750.15 


2.42 


Fig.|A3| 


Nil 


6767.76 


1.82 


Fig.lAH 



(1): The "standard" bisector line. (2): Strong blending evident in bisector, not plotted. 



there are aspects of the bisectors in Fig.^that call for attention, 
in particular the peculiar appearance of the y Ser bisector, and 
the shape of the EK Eri bisector. 

As can be seen, the Fe 1 6252 line is quite weak in y Ser 
Moreover, the lines are broadened by rotation, which further 
reduces the depth of the line, but on the other hand allows a 
better sampling of the profile. However, the star is close to the 
granulation boundary so no definite conclusions can be made. 

EK Eri is a well-known over-active but slowly rotating star, 
exhibiting long-period photometric and spectr oscopic varia- 
tions due to rotational niodulation of stars pots ( IStepieniri993t 
IStrassmeieret al.!' 19991: ball et aljE005al) . The bisector does 
not much resem ble the "canonical" shape of a G8 giant (e.g. 
rj Dra from lGravl2005[ Fig. 3), but rather indicates a somewhat 
hotter spectral type, especially towards the line core. Since the 
core is formed in the higher atmospheric layers, this could in- 
dicate a substantial heating of the upper layers by the chro- 
mospheric activity. Fitting the bisector to the tabulated ones of 
Gray, we find that the lower part is better described by a G5 
star of luminosity class III or II (no perfect fit can be obtained), 
indicating a heating of the upper layers of the photosphere by 
~ 200-900 K. This is however quite speculative. Since the star 
is known to have spots, the bisector shape may be reflecting the 
contributions of different velocity fields in the spots. It could 
also simply have a misclassified spectral type. 

We will now proceed to an investigation of other potentially 
useful isolated spectral lines (Table|3. These lines were chosen 
based on VALD linelists (Kupka et al. 1999) of the solar spec- 
trum', by requiring that they be unblended, isolated and rea- 
sonably strong. The bisectors of the Car 6499 line were found 
to be very irregular for all our sample stars, although no ob- 
vious blends could be identified, and they are not plotted. The 
bisectors of the remaining lines are shown in the Figures in 
AppendixlAl 

The Fez /i 615 1.62 line is very well defined, but not very 
deep in all our stars. The line may be influenced by weak con- 
tributions of Nil /16151.67. The closest major line however, is 
more than 0.5 A away (Si i /16152.29). Examining Fig. lA.ll we 
find no obvious systematics, and we consider this line inappro- 
priate for stand-alone analysis. This seems to indicate, that the 
presence of even very weak blends can disturb the interpreta- 
tion of the bisector, provided the lines are sufficiently close in 
wavelength. 



' The Vienna Atomic Line 

|http : //'ams . a str o . univie . a c . at/vald/ 



Database, 



The Fei /16252.56 is the standard bisector line, thought to 
be sufficiently well isolated and free of blends. For low rota- 
tional broadenings, the neighboring Vi /16251.83 line is well 
separated, and the only other weak lines that may possibly con- 
tribute are all on the order of 0.5 A away or further, namely 
Rui ^6252.07, and Cri /16253.17. 

As mentioned, the bisectors of Cai /16499.65 were all 
very irregular This may be ascribed to blending with nearby 
Si I ^6499.23 and N i 6499.5 1 . 

The Nil /1 6643.63 line is very deep and well defined for 
most of our stars, and generally shows a larger velocity span 
than the Fez 6252 line. It is probably affected by blending 
with Vi /16643.79 and Gdz /16643.94 on the red side and with 
Srz /16643.53 and Crz /J6643.23 on the blue side, which may 
contribute to the larger velocity span, but rendering any physi- 
cal interpretations very uncertain. 

For the Fez /1 6750. 15 line, only a few of our stars show 
bisectors that resemble the classical "C" shape (Fig. IA.3> . 
The line and bisector are well defined, but the sharp redward 
turn in the upper part may indicate blending with another 
line. The prime candidates for blending are Fez /16750.58 and 
Cu z /i6749.65, although the latter is probably too far away. 

The Niz /1 6767.76 line show well defined bisectors for 
most stars (Fig. lA.4l l. However, no obvious systematics in the 
shape makes it unlikely that this line can be used to infer atmo- 
spheric parameters. It is probably disturbed by Fez /i6767.71, 
and possibly by Coz /16768.17. 

In conclusion, it seems that a "well-isolated line" translates 
to a minimum wavelength separation to potential blends of ~ 
0.5 A. 

5. CCF bisectors 

It was shown by iGravl ll2005h that the blue-most point of the 
classical single-line bisector is a very good indicator of lumi- 
nosity class for the late G to early K stars. If an analogous rela- 
tionship exists for the CCF bisector, then it will provide us with 
a luminosity estimate for stars where high-S/N spectroscopy is 
not feasible. 

The mean CCF bisectors are plotted in Fig. |2l Their ap- 
pearances differ somewhat from the Fe z 6252 bisectors, which 
is not surprising given that they are the bisectors of a "mean" 
spectral line, and given the diversity of individual line bisectors 
evident from Figs. lA.llA.4l 

It must be stressed that one must be careful in comparing 
CCF bisectors of stars with different spectral types, since the 
bisector will depend heavily on the mask used in computing 
the CCF. An example of this is shown in Fig. |3lfor EK Eri, 
using G2 and K5 masks. Hence, one would suspect that in 
general the CCF bisectors can only be used to compare stars 
of similar spectral types, with CCFs computed using the same 
mask. However, as we use a G2 mask to calculate the CCF for 
all stars except a Hor, this precaution is of less importance in 
the present study. Another concern about the masks, is whether 
they are constructed in a way that optimises the physical inter- 
pretation of the CCF bisector These masks are constructed to 
yield the best possible RV precision, which means excluding 
very weak, very strong, and heavily blended lines. So while a 
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EK Er 
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J^ 
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Q M e r 
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Fig. 2. The mean CCF bisectors of all program stars, grouped 
according to luminosity class. The vertical dashed lines indi- 
cate the RV of each star. Formal eiTor bars are the size of the 
plotting symbols or smaller. 



E EK Erl (GSIII) 
0,9^ 




6.95 7.00 7.05 7.10 7.15 

RV [km/s] 

Fig. 3. The CCF bisector of EK Eri computed using two differ- 
ent masks. 



given mask will include lines of many different elements, it will 
be free from strong blends and hence represent a near-optimal 
"average" spectral line. 

Returning to Fig.|2 there seem to be trends in the bisector 
shapes, possibly with temperature, but most noticeably with lu- 
minosity and/or gravity. In Fig.0]we plot the absolute magni- 
tude as a function of BIS + Cb, which seem to exhibit a well- 
defined linear relationship. A similar relationship is found plot- 
ting log^ versus BIS + ci,. We derive the approximate relations 



M, = 3.86 - 24.88(BIS + ct) 
logg = 4.06-9.22(BIS+Ci), 



(3) 
(4) 



which we must stress are only approximate, are based on very 
few data points, and should not be used to derive final stellar 
parameters. 

If proven to be correct, these relationships can be used to 
measure the luminosity from the shape of the average CCF bi- 




BIS + Cb [m/s] 

Fig. 4. Absolute magnitude as a function of BIS + Cb- 

sector, in analogy with the blue-most point luminosity indicator 
for the classical bisector, only that these relationships may be 
extendable to all the late spectral types. However, to establish a 
reliable relationship, or indeed to determine whether such a re- 
lationship exists at all, more data is needed. Obviously, the fact 
that the values of the bisector measures depends on the mask 
used, calls for extreme caution. A more thorough analysis of a 
larger sample is cuiTently in progress (Dall et al., in prepara- 
tion). 

6. Variations of CCF bisectors 

In the remainder of this paper we will concentrate on the 
changes in the CCF bisector measures, rather than the absolute 
values of these, and examine how these changes may be used 
as diagnostic tools when looking at only one star. In doing so, 
the dependence on the mask used for the CCF disappears, as 
we look only at relative changes. 

First, however, we must address the expected scatter of the 
measures in the case where no activity or oscillations are likely 
to show up, using the non-active stars with low-amplitude oscil- 
lations, i.e. primarily yU Ara and the Sun. Some of the remaining 
targets will then be addressed in turn. 

A general idea of the expected scatter can be gained by 
looking at Fig.|5] which shows the scatter in the bisector mea- 
sures for all the stars. Obviously, the noise will depend on the 
S/N of the CCF and on the spectral type. As can be seen, ex- 
cept for y Ser and EK Eri, the scatter is very similar over the 
full range of spectral and luminosity classes. In particular, the 
RMS spread and the point-to-point scatter are almost identical. 

6.1. No correlations -/i Ara 

H Ara is the host of a low-mass planet orbiting the parent star 
in 9.5 days dSantos et alJl2004l) . Moreover, it exhibits well- 
described low-amplitude solar-like p-mode osc illations with 
individual mode amplitudes up to 0.4 ms^' ( iBouchv et alJ 

I II 1 ^ 1 r 

2005; Bazot et al. 2005). From the ESO Science Archive we re- 
trieved 100 consecutive CCF observations, spanning 3.7 hours 
— long enough to include several oscillation cycles and short 
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Fig. 5. The scatter in all four bisector measures for all program Fig. 7. The Sun: The bisector measure BIS as function of the 



stars. Note that y Ser is close to the granulation boundary, and 
EK Eri is an active star (discussed in the text). The other sym- 
bols are for t Cet (O), Sun (A), /3 Aql (+), jj. Ara (n), (3 Hyi (x), 
6 Pav (•), a Men (>) and a Hor (o). The individual Icr RMS 
scatter (upper point) and point-to-point scatter (lower point) are 
connected with dotted lines. 



angular distance to the sun 9 for one 2.6 hour segment. The 
same trend is observed in bb, while no such trends can be seen 
for Vbot and Cb- 



6.2. Scattering of sunlight in the Earth atmosphere 
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Fig. 6. The bisector measures as functions of RV for// Ara. The 
scatter is dominated by solar-like ;?-mode oscillations; see text. 
Horizontal axis is RV in all plots, with the scale indicated in 
the upper right panel. 



The spectral lines in solar spectra obtained on the d aytime sky 
-^° are kn own to exhibit EW and line depth changes. lOrav et alJ 
_ (yOOO) measured the line depth variations and found them to 
^ be symmetrical, i.e, they found no changes in the bisector due 
g(j> to asymmetries of the lines, as a function of the angular sep- 
^ aration 9 between the telescope pointing and the position of 
-70 ^ the Sun. Their detection limit on the bisector variations was 
20 m s ' , and they mention the possibility of smaller variations 
7 being present. They successfully modeled the line depth varia- 
tions in terms of aerosol and Rayleigh-Brillouin scattering as a 
function of 0. 
^ < Our data, taken in directions close to opposite the sun, cov- 
6 7 ers a span in of ~ 30°, in which we observe a smooth change 
^ in the CCF line depth on the order of 1 %, consistent with the 
findings of lOrav et alJ yOOO). From our data it is evident that 
there is an asymmetric contribution to the atmospheric scatter- 
ing, as revealed in the BIS, which varies with 4 ms"' over the 
interval (Fig.0. Obviously, a better sampling of from 0° to 
180° would be desirable, as we demonstrate that line asymme- 
tries in the scattering can indeed be detected with HARPS. 



enough to avoid longer-term effects from the presence of the 
planet. 

The four investigated bisector measures are shown in Fig.|6l 
The ~ 10 ms ' scatter in RV (horizontal axis) is dominated 
by the p-mode oscillations. Note that the measured velocity 
is a sum of a few dozen modes, each of which has small 
(~ 0.4 ms"') amplitude. These modes beat to produce t he ob- 
serve d peak-to-peak velocity excursions (see e.g. Bouch v et al.l 
I2OO5 . Fig. 5). It is evident that no obvious correlations are 
present, and we will then adopt the scatter shown in these plots 
to be "typical" of non-correlating bisector measures. 



6.3. Variations due to oscillations? 

The stars a Men, 6 Pav, /3 Hyi and a Hor are all ex- 
pected to show solar-like p-mode acoustic oscillations (see 
iKieldsen et alEo OS. for 6 Pav and /? Hyi). 

Our data on a Men were taken in individual six-spectra seg- 
ments of less than one hour duration, over a period of three 
months. Peak-to-peak variation is ~ 6 ms ', and no correla- 
tion with bisector measures can be seen. 

The observations of 6 Pav span 1.8 hours, and show 
clear solar-like oscillations with peak-to-peak amplitudes of 
~ 10 ms"'. No correlation with bisector measures can be seen. 
Higher scatter of the bisector measures compared with the other 
stars is evident from Fig.|5l We have no explanation of this fact 
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other than the bad weather conditions (bad seeing and wind 
from the south) prevailing at the time of observation. This, and 
the short exposure times, could also be part of the reason for - 
the low S/N in the combined spectrum. 

The results of the asteroseismic campaign on p Hyi will be ' 
published elsewhere (Kjeldsen et al., in preparation). Here it 
suffice to note the peak-to-peak amplitude of the /3 Hyi oscil- 
lations of up to ~ 10 ms"', revealing weak coiTelations with 
bisector measures BIS and b/,. A Fourier analysis of the bisec- 
tor measures reveals that the oscillation signal is clearly present^ 
in the BIS and bb data, while Vbot and ci, do not reveal any oscil- 1 
lation signal at the 6cr level. This could indicate that the upper 
part of the line is insensitive to the oscillations. 

a Hor is one of the most stabl e Hipparcos sta rs, with a 
photometric error of < 0.3 mmag (Adelman 2001). No reli- 
able stellar atmospheric parameters c ould be found in the liter- 
ature, so we performed an ATLAS9 ( lKurucall993l) model at- 
mosphere fitting using Fe i/Fe ii ion balance, following the pro- 
cedure of Dall et al. (2005b) . The RV curve of a Hor shows 
a smooth variation over the duration (~ 2.5 h) of the obser- 
vations. The peak-to-peak amplitude is ~ 25 ms"', maybe ^ 
slightly larger as we have not covered the full period. Assuming 
this to be due to stellar oscill ations, using the scaling relation of 
iKjeldsen & Bedding I (119951) predicts a photometric amplitude 
of ~ 0.3 mmag, i.e. the oscillations should have been barely 
detectable by Hipparcos. More importantly for our purpose, 
variations at this level cannot be seen in the bisector. 

In summary, solar-like p-mode oscillations have a weak in- 
fluence on the bisector, and are only properly revealed by a 
Fourier analysis. 

6.4. Variations due to activity 

The star EK Eri shows a clear example of activity induced 
bisector variations. The activity and associated variability of 
EK Eri has been studied extensively for years, and it is gener- 
ally agreed to be a slow rotator with a fossil Ap-type magnetic 
field (Stepieii 1993; Strassmeier et al. 1999). Dafl et al. (2005ai 
found the BIS of the CCF to correlate well with the RV and the 
activity index log R[i^, meaning that the activity displays a co- 
herent pattern from the photosphere to the chromosphere. 

We examined all four bisector measures, shown in Fig.|S] 
in order to test which is best suited to describe activity induced 
RV changes. In addition to the BIS, bt and Vbot show clear cor- 
relations with the RV, while c/, seem to be insensitive to the 
activity, which is also evident from Fig.|5] Hence, it seems that 
the stability of c/, is real, and that the curvature of the CCF 
bisector is unaffected by activity. Note that the slopes of the 
correlations are positive for BIS and bi, as expected from activ- 
ity and negative for Vbot, and that the peak-to-peak amplitude is 
of the order ~ 100 m s' , i.e., much higher than the oscillation 
scatter. 



6.5. Variations due to unseen companion? 

For y Ser, the scatter in all the bisector measures are high 
(Fig-Ell. A. closer inspection reveals clear correlations with BIS 
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Fig. 8. The bisector measures for EK Eri. Horizontal axis is RV 
in all plots, with the scale indicated in the upper right panel. 




Fig. 9. The bisector measures as functions of RV for y Ser. 
Horizontal axis is RV in all plots, with the scale indicated in 
the upper right panel. 



and bb, less clear with Vbot and no correlation with Ch (Fig.|5}. 
Note that the slopes of the correlations are opposite to the case 
of the activity induced variation of EK Eri. This hints that the 
variation may be due to either a fainter orbiting star or a chance- 
alignment with a background star, corresponding to the case of 
HD 41004 (.Santos et aL.2001) where the blending with light 
from a much fainter companion produced a negative slope of 
BIS versus RV. Such a coiTelation could also result from chance 
alignment with a fainter background star, modulated by seeing 
variations. Note that the possible companion discussed here is 
a stellar companion, contributing flux to the total spectrum. A 
planetary companion would not show up in this way. 

We hesitate to propose that the correlations could be due 
directly to oscillations, since the apparent amplitude of the RV 
variation is lower than for a Hor, where no correlations were 
found, and similar to /3 Hyi, where only weak correlations were 
found. However, the slope of the correlations are the same as 
for /3 Hyi, so this may actually reflect oscillations, and the rea- 
son for it showing up in y Ser and not in a Hor may lie in the 
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proximity of the former to the granulation boundary. This ques- 
tion must be resolved by observing a larger sample of stars. 



7. Concluding remarks 

In this paper we have examined a small sample of well-studied 
bright stars of different spectral types and luminosity classes. 
We have measured the bisectors of isolated spectral lines using 
high S/N combined spectra, and we have measured bisectors of 
CCFs of individual exposures, obtaining a high S/N from the 
combination of hundreds of individual lines. 

The main purpose of this paper — to relate the single-line 
classical bisector to the single-spectrum combined-line CCF 
bisector — has been achieved: We have shown that one may 
use the CCF bisector in much the same way as one would use 
the single-line bisector. Moreover, we have shown that it is pos- 
sible to use the defined CCF bisector measures for quantitative 
analysis, employing both the absolute values and the variations. 

For the single-line classical bisectors, we can point out the 
following: 

- We confirm the general shapes of the bisectors of the 
Fe 1 6252 line for stars of di fferen t spectr al and luminos- 
ity classes, as published by lOravl (120051) . using sums of 
HARPS spectra to achieve high S/N. We did not find other 
spectral lines with similar well-behaved bisector patterns 
with spectral type. We find that to be well isolated, a line 
must not have any potential blends closer than ~ 0.5 A, 
based on the six lines we have investigated. 

- Based on the bisector shape of the Fe 1 6252 line in the spec- 
trum of EK Eri, we suggest that this may point to an ex- 
tra heating of the upper photosphere of several hundreds 
degrees, presumably supplied by the active chromosphere. 
However, the peculiar bisector shape may also be due to 
photospheric spots, or to inaccurate spectral classification. 

For the CCF bisectors, we can point out the following: 

- We find that the CCF bisector could be used to approxi- 
mate the luminosity and log g using the absolute values of 
the bisector measures BIS, bh, v'bot, and c/,, and we derive 
approximate relationships. It is slightly surprising that a lu- 
minosity relationship can be recovered in the CCF bisector, 
since it includes many different spectral lines with different 
individual bisector shapes, not only between lines of dif- 
ferent elements, but between lines of the same element too 
(Fig.^and the Appendix). An added advantage of the CCF 
bisector in this case, is that the relationship seems to apply 
to all the spectral types of our sample, and is not dependent 
on the existence of a blue-most point, which disappears for 
the hotter stars. Naturally, the validity and accuracy of such 
relationships will be the subject of future study of a larger 
sample of stars. 

- We find that normal p-mode solar-like oscillations do not 
reveal themselves in the bisectors at a significant level (for 
peak-to-peak RV amplitudes < 25 ms '). However, there 
may be a dependence of the bisector response to oscilla- 
tions on the spectral type or luminosity class. For y Ser 



we do find correlations between RV and the bisector mea- 
sures BIS and bh, even though the oscillations are at a very 
low level. We speculate that this could be due to an unseen 
stellar companion, either a physical companion or a chance 
line-of-sight alignment. Another possibility is that y Ser is 
too close to the granulation boundary to allow any mean- 
ingful interpretation of its bisectors. 
- For EK Eri, the only active star of our sample, we find the 
RV correlated with all bisector measures except the curva- 
ture Ch- Note that the slopes of the correlations for y Ser 
and EK Eri are reversed, giving an indication of different 
mechanisms being involved. 

Finally, we confirm the effects of aerosol and Rayleigh- 
Brillouin scattering on solar spectra taken on the daylight sky. 
Moreover, we show for the first time that the bisectors and 
hence the lines change shape, thus there is an asymmetric com- 
ponent of the scattering. 
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Appendix A: Bisectors of additional lines 
List of Objects 



51 Peg! on page 2 
Sun[_on page 3 

FMf on page 3 
^Cef on page 3 
^Hoifon page 3 
E^^rif on page 3 
^Men| on page 3 

Sen on page 3 
u Aral on page 3 



Am on page 3 

Paa on page 3 

77 Drai on p age 4 



np; 



HP 41004 on page 7 



